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Similarities between models of fragmenting nuclei and disordered systems in condensed matter 
suggest corresponding methods. Several theoretical models of fragmentation investigated in this 
fashion show marked differences, indicating possible new methods for distinguishing models using 
yield data. Applying nuclear methods to disordered systems also yields interesting results. 
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I. INTRODUCTION 

The breakup of nuclei into clusters of various sizes in 
high energy heavy-ion collisions is currently being studied 
both experimentally and theoretically |§-§< The 

reasons for undertaking such studies are many, but one 
is to determine whether a liquid-gas like phase tran- 
sition occurs at densities and temperatures away from 
the normal conditions of the nucleus (achieved, albeit 
briefly, during the collision). By investigating the behav- 
ior of particular functions of the fragmentation pattern, 
the character of the transition can be determined, and 
the various critical point exponents can be measured. 
For example, the EOS collaboration J^|| studied a set 
of moments introduced by Campi pOf to determine the 
three critical exponents /?, 7, r and the critical multiplic- 
ity. Their results strongly suggest that nuclear fragmen- 
tation can be understood in terms of a liquid-gas or pos- 
sibly a percolative phase transition. 

Many non-nuclear systems also exhibit clustering. For 
example, in condensed matter physics the phase space 
of disordered systems [ pd]JT2| such as spin glasses [ ^3|Jl4| 
exhibits clustering of states of the system around energy 
minima. As in the nuclear case, the sizes and distri- 
butions of the these clusters are studied in an effort to 
better understand the character of these systems and the 
nature of the disorder. As such one is not surprised to 
discover a great degree of overlap in studies of nuclear 
fragmentation and disordered systems. Differences in the 
approaches are considerable however. 

This paper investigates in detail the correspondences 
between fragmenting nuclear and disordered systems, dis- 
tinguishing the similarities and the differences in the 
approaches. A similar investigation by Higgs |i~5| ] dis- 
cussed such parallels between biological and disordered 
systems. Our purpose is to enlarge the avenues of in- 
vestigation in nuclear fragmentation by examining areas 
where parallel methods of analysis are likely to be fruit- 
ful. Such approaches have been successful before, such 
as the application of percolation theory to nuclear frag- 
mentation p6|JTc| ] and our preliminary results presented 
here and elsewhere |l7j encourages such a detailed look. 
Our principal result, that all the models considered can 



be realized as a sequential breaking of the interval, rec- 
onciles the observation of these similarities to a distinct 
mechanism. 



II. FRAGMENTATION OBSERVABLES 

A nucleus of A nucleons (or in general a system of 
A objects) can be partitioned into n k clusters of size k. 
The total number of clusters (in nuclear physics known 
as the multiplicity) then is given by to = and the 

total number of nucleons A = ^ fe kn^. The set of n^'s 
subject to knk — A determines the partitions of A, 
with n = (ni, n%, ■ ■ ■ tia) describing a particular partition 
or fragmentation. Alternately, one can describe such a 
partition by A = (Aj., . . . , A m ) where A& is the size of the 
&;th cluster, A = A&, and the clusters are ordered in 
some fashion, e.g. according to size Ai > A2 > ■ ■ • > A m 
or for dynamical fragmentation by order of appearance. 

Either n or A describes the clustering pattern, but we 
are more interested in functions of the pattern which elu- 
cidates the character of the fragmentation. The moments 
of a fragmentation pattern are often examined for this 
reason, as they distill the results of a complex fragmen- 
tation pattern into a single number which can be simply 
interpreted. We define the sth moment of a fragmenta- 
tion pattern n by 



Y s (n) = yk s n k 



k 



(1) 



The first two moments count the number of clusters and 
objects, i.e. Yq = m, Y% = A. The second moment 
Y-2 = Y indicates the average size of the cluster (Y/A) a 
randomly chosen object belongs to, and higher moments 
can be similarly understood. In spin-glasses, Y/A 2 is 
known as the participation ratio in Anderson localiza- 
tion Q. 

Sometimes we wish to reduce these moments Y s , elimi- 
nating any contribution to them from the largest cluster, 
i.e. we define the reduced sth moment as 



M,(n) = r s (n)-fc^ 



(2) 
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where fc max is the size of the largest cluster. Such a 
moment seems natural if we consider the system as a 
liquid-gas system. Then the largest cluster can be con- 
sidered the liquid phase, and such a reduction corre- 
sponds to considering the contributions from the gas va- 
por only. Alternately, percolation theory would identify 
such a cluster with the incipient infinite cluster, and the 
reduction would separate the system into percolating and 
non-percolating regions. In spin-glasses, such a reduction 
removes the contribution of the extended state, the deep- 
est valley in the rugged energy landscape. 

III. FRAGMENTATION MODELS 

The discussion of observables considered only a sin- 
gle fragmentation event. Physical systems produce such 
events in great numbers, but not necessarily with a uni- 
form distribution. For the observables discussed above, 
once the distribution of events is specified any ques- 
tion about the observable can in principle be answered. 
Clearly there is a physical mechanism which leads to such 
a distribution, but it is possible that several mechanisms 
may yield the same fragment distribution. Analyzing en- 
semble averages or distributions on fragment observables 
would then fail to distinguish between various mecha- 
nisms. Therefore, one can specify a model by its un- 
derlying fragment distribution without recourse to the 
mechanism which generates it. This section does this for 
a number of cases applicable to nuclear fragmentation 
and disordered systems, and shows that all these systems 
can be generated by a random breaking mechanism. 

A. Nuclear Systems 

When two heavy ions collide at intermediate energies, 
the nuclei dissolve into a diverse mix of fragments. The 
fragment distribution is necessarily complex, but its gen- 
eral statistical features are determined by simple thermo- 
dynamics. For instance, the kinetic motion of each frag- 
ment limits the available phase space, and so the proba- 
bility of seeing m fragments is proportional to x m , where 
x = V/Xj,, the ratio of the available volume to the ther- 
mal volume. Internal excitations and other effects will 
also enter, but in the large T limit kinetic motion domi- 
nates, predicting a large numbers of fragments. Low tem- 
perature considerations likewise constrains our fragment 
distribution. At low T a large fragment accompanied by 
some very small fragments is expected. The region be- 
tween the two natural limits having several intermediate 
mass fragments is more difficult to describe and is the 
subject of current research. 

One might believe that the distribution of fragments 
in the intermediate region is not dominated by thermal 
factors, but rather by the available phase space associ- 
ated with the number of possible patterns. Sobotka and 



Moretto [[L9| explored the simplest such model in which 
every pattern is given equal probability. Such a model 
shows an exponential fall-off in the fragment yield, con- 
trary to the experimentally observed power law. One can 
generalize this model in many ways, the most obvious of 
which is to consider a model where different cluster sizes 
would contribute different amounts to the fragment dis- 
tribution phase space, i.e. Pr(n) = (1/Za) 11/,- x k k so 
that Xk — 1 reduces to Sobotka and Moretto's model. 
The partition function Za = Snllfc^fc'" m this case 
can be calculated using the recursion relation Za = 
( l / A )J2k>a k J2i>o x k Z A^k with Z = 1. 

Such models have many interesting properties, but 
need some modifications to avoid problems associated 
with the indistinguishability of clusters of the same size. 
This indistinguishability should reduce the size of the 
phase space available to the fragments, or equivalently, 
reduce the probability for a particular fragment distri- 
bution by factors of l/n^l after the argument of Gibbs. 
With this consideration, a simple and somewhat generic 
model of fragmentation would be to take the probability 
of a particular fragmentation pattern to be given by 

Pr(n) = J- J] ^ (3) 

Za fc>o nkl 

where the partition function or normalizing factor Za can 
be obtained using the recursion Za — (1/^4) Ylk kxkZA-k 
with Zq = 1. Such a model as applicable to nuclear 
fragmentation was explored in | j20|]2lf| , and took Xk = 
xy k ~ 1 J '0k where x = V/X^ arises from the thermal mo- 
tion, y describes effects due to binding and internal en- 
ergies and (3k allows different sized clusters to contribute 
different amounts to the phase space. Thermodynamic 
arguments discussed in ]22| allow y to be expressed as 
y = exp{a v /k B T + (kBT/e )(T Q /(T + T ))} where a v 
is the binding energy per nucleon (i.e. it is assumed that 
Esik) ~ av{k — 1) for a cluster of size fc), e is the 
energy spacing between excited levels and To is a cut- 
off temperature for internal excitations. The parameters 
(3k can be taken as k T so as to reproduce the observed 
cluster distribution, (nk) ~ fc~ r , where r is a critical 
exponent introduced by Fisher |23]| . For nuclear frag- 
mentation r w 2.2 Q. 

Several models with (3k — k T have been explored by us 
and others. For example, (3k = I, called the linear chain 
model, was extensively studied by Gross et al. p6[ . 
Another example, (3k — fc, was the first model investi- 
gated by one of us (A.Z.M.) p7|] . We have investigated 
(3k = k T with t > 2, in particular t = 5/2 more re- 
cently |20| ]. It has many features in common with Bose- 
Einstein condensation and percolation theory; namely 
below a certain value of cluster forms con- 

taining a finite fraction of the mass of the system even in 
the infinite A limit. This "infinite" cluster can be likened 
to the accumulation of particles in the ground state in a 
Bose gas below T — T c , or to the formation of the incip- 
ient infinite cluster in percolating systems above p = p c . 
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Related weights including more factors (some dynam- 
ical in nature) have been studied by Gross et al. |6| and 
Bondorf et al. H in describing more detailed models of 
nuclear fragmentation. 



B. Disordered Systems 

Disordered systems can also involve clustering, and in 
those cases the probability of seeing a particular cluster- 
ing pattern is of interest. In the case of nuclear fragmen- 
tation, thermodynamic and phase space arguments were 
used to arrive at the pattern distribution, i.e. Eq. (||). 
Here we arrive at a similar description by analyzing a sim- 
ple but general mechanism which generates disorder: the 
random sequential breaking of an interval. This model 
was studied extensively by Derrida and Flyvbjerg [ p8|j29| 
and more recently by Frontera et al. fl30| |. Here we gen- 
eralize their mechanism in such a way as to encompass 
Eq. (Q) and many other models. 

The sequential breaking of an interval is simple to de- 
scribe. In the continuous case, the unit interval is se- 
quentially broken into pieces of sizes W%,W2, ■ ■ ., each 
chosen randomly in some fashion. Specifically, the se- 
quence Wi = zi, W 2 = (l-z 1 )z 2 , W 3 = {1-Zi)(l-Z2)z 3 
et cetera, is generated where each Zk is chosen from some 
probability distribution pk(z). In other words, at step i 
one breaks off a randomly chosen fraction Zi of the re- 
maining piece of size 1 — W\ — . . . Wi-\. Typically the 
probability distribution on Zi is independent of i, though 
that need not be the case. Derrida and Flyvbjerg studied 
the case where p(z) = x(l — z) x ~ x f2S|j , 

Let us consider an equivalent, but discrete model. In 
this model, an interval of length A is partitioned sequen- 
tially into integral sized pieces. Suppose a piece of size k 
is broken off from the interval of size A with probability 

where (x) n — x(x + 1) • • • (x + n — 1) and x > is a free 
parameter. In the limit A — > oo, k/A — > W, this distribu- 
tion converges to the continuous distribution mentioned 
above. 

Applying this process repeatedly until the interval is 
exhausted, the sequence Ai, A2, . . . , X m is generated. The 
overall probability for this sequence appearing from this 
process is the product Pr(Ai; ^4) Pr(A2; A — Ai) Pr(A3; A — 
Ai - A 2 ) • • • Pr(A m ; A - X± - ... - A m _i) and is given by 



particular unordered clustering pattern appears is deter- 
mined, and is given by 



Pr(Ai,...,A m ) =x r - 



A + x - 1 
A 

1 



A(A - Ai) ■■ - (A — Xi 



— A r , 



1) 



(5) 



Pr(n) 



A + x - 1 
A 



n ?lfc i (k) 



(6) 



By summing over all possible orders of such a sequence, 
the probability Pr(n) = Eamii P r (^iv)^m) that a 



where we have applied Eq. ( Al). This is in fact the distri- 
bution given by Eq. (^) with Xk = x/k. More generally, 
the sequential breaking process where the probability of 
breaking off a fragment of size k from an interval of size A 
is given by Pr(fe; A) = (kx k /A)Z A - k /Z A yields Eq. (|) as 
the fragmentation pattern. This gives a possible mecha- 
nism for obtaining such a fragment distribution, but it is 
not unique. Another mechanism based on Markov chains 
was explored by us before 

By modifying the random breaking mechanism we can 
obtain the other model discussed in the previous section 
(the equipartition model) as well. Suppose one allows 
multiple clusters of the same size to be broken off in- 
stead of just one. For example, take the probability that 
to clusters of size k are broken from an interval of length 
n as Pr(k,m;A) = (kxf /A)(Z A - m k/Z A ). Then the dis- 
tribution of cluster sizes can be shown to be given by 
Pr(n) = (^/Z A )Ylk x7 k k ■ This is in fact a generalization 
of a well-known method for generating equiprobable par- 
titions introduced by Nijenhuis and Wilf ]3~i| . 

Having obtained the models of the previous section as 
cases of this disordered system mechanism, we are en- 
couraged to believe that other disordered systems can 
be likewise described in the context of randomly break- 
ing the interval. Indeed, two other models considered by 
Derrida and Flyvbjerg fit nicely into this mechanism. 

Consider the quenched random map model J2^,|32|,|33| . 
The A A maps or functions of A points onto themselves 
are chosen uniformly at random, and their cluster struc- 
ture is determined by the number of points in each basin 
of attraction, where a basin is defined as a limit cycle of 
the iterated action of the map. For example, the map 
1 h 2,2 h 1,3 h 2,4 h 3 has a single limit cycle 
(1 1— » 2 1— ► 1 . . .) which all the points eventually enter, 
so that the cluster structure is just one single cluster 
of size four. Such a simple model has interesting prop- 
erties, and its cluster pattern is given by Eq. (|) with 
Xk = Sj=o j ]\ ~ e k /2k. This can be inferred from 
the fact that the probability that the map is indecom- 
posable is known and would be given by Xk/Zk where 
Zk = k k jk\ by simple counting. Interpreting Eq. (|j|) in 
terms of randomly breaking an interval has already been 
accomplished, so we see that the random map model is 
a special case of randomly breaking the interval. The 
asymptotic limit can be seen by considering the elements 
in the sum as terms in the Poisson distribution, e~ x x J / ' j\ 
with x = k. The factor e k does not affect the cluster dis- 
tributions, so that the model asymptotically converges 
to Xk = l/(2fc), and hence is included in the original 
distribution considered by Derrida and Flyvbjerg |2§| ] 
p(z) = x(l - z)*- 1 with x = 1/2. 
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Secondly, consider what happens when the sequential 
pieces are chosen from a changing distribution. For ex- 
ample, suppose the i+lst piece is chosen from the interval 
of size Ai = A — Ai — ... — Ai, Ao = A with probability 

\k-l) (x + 7 + i(l - 7))a 4 -i V ; 

In this case the sequence of pieces is generated with prob- 
ability 



Pr(A 



1 3 ... 3 A m 



Xm) = A\ 



{ x ) (m-l;l--y) 

(x + j)a-i 



1 



A(A - Ai)(A - Ai - A 2 ) • • ■ A fc A (A*, 

where X( n . Q ) is a generalized Pochhammcr symbol 

__ J 1 for n = 

»(n;a) = j a;( x + a ) . . . ( x + ( n _ 1 ) a ) for n > 

Summing over all the possible orders of the A's, gives the 
probability for the unordered fragmentation pattern 



Pr(n) = AI 



(a0(m-l;l-7) 

(x + j)a-i 



n nfe ! 



(7)fe- 
fc! 



(9) 



This is clear since the only term that depends on the 
order of the fragments is \/{A(A — A\) ■ ■ •); the other 
terms are invariant under perm uting the order of the A's, 
so that we can apply Eq. flAl| ) directly. 

This distribution pattern has been studied in a statis- 
tics context by J. Pitman p4| , but what we would like to 
emphasize here is that in the limit x + 7 = 1, this model 
reproduces asymptotically all of the clustering features 
of a Sherrington-Kirkpatrick spin-glass p5| , ^6[p9| . Un- 
like the other disordered systems, we cannot prove this 
result, but the evidence (discussed below) is compelling 
enough that we consider Eq. ([)[) with x+j = 1 as a model 
of a spin-glass. Why the changing distributions produces 
the patterns seen in a spin-glass is perhaps unexpected 
and should be explored further. 



IV. TRADITIONAL ANALYSIS OF MODELS 

Now that we have a description of the fragment dis- 
tributions we are in a position to make definitive pre- 
dictions. Since nuclear and disordered systems have 
spawned separate methods of analysis, we begin by 
briefly considering some features of these models from 
their traditional points of view. 



A. Nuclear Systems 

Fragmenting nuclear systems are traditionally ana- 
lyzed by studying expectation values of the various ob- 
servables as functions of thermodynamic variables. Some 



of the important observables are the multiplicity m and 
its variance, the size of the largest cluster fc max , and the 
reduced moment M 2 . They are important because of 
their expected behavior in the region of a critical point 
phase transition. For example, from percolation theory 
we know that the reduced moments are functions of the 
critical exponents and p — p c in the region of the critical 
point p|, i.e. 



(M ) = (m)-l~\p-p c \ 2 - a 
(M 1 )=A-(k nm ) ~ \p-Pcf 
(M 2 ) = (Y 2 )-(k 2 max )~\p-p c 



(10) 



The EOS collaboration || has argued that both nuclear 
fragmentation and percolation theory can be analyzed 
in a similar fashion, extending the observations and re- 
sults of Campi 10| 37]] . Thus from observing the behavior 



of the moments the critical point and exponents can be 
deduced, and a number of papers in the literature have 
studied this in a variety of theoretical |37]|38|,^6| and ex- 
perimental systems. 

In the model given by Eq. (|^), these quantities can 
be computed from their definitions and by utilizing the 
following results 



(nit) = x k 



ZA-k 



(nj(n k - 5jk)) = x 3 x k 



ZA-j — k 

Za 



(11) 

(12) 



Computation of the reduced moments additionally re- 
quires that various expectation values of fcmax be deter- 
mined. This is most easily done by expressing the expec- 
tation values in terms of the probability distribution on 
fcmax- The probability that a random fragmentation of 
nucleons has a cluster of size fc ma x &s its largest cluster is 
given by Pr(fc max ) = AZ A (x 1) x kmBX )/Z A {xi, ...,x n ) 
where AZ A (xi, . . . , #fc ma J = Z A {x\, . . . , x femax , 0, . . .) - 
Z A (xi,. . . ,x femax _i,0, . . .). When fc ma x > A/2 a simpler 
expression can be used. Since any cluster of size k > A/2 
is automatically the largest one can show that AZ A = 
x k Z A - k , so that Pr(fc max = k) = x k Z A - k (x)/Z A (x) 
when fc > A/2. 

For typical models such as x k = x/k T with t > 2, the 
variance of m and the reduced moment M 2 should peak 
at the point x c /A = 1/£(t— 1). The finite size effects can 
be considerable, and the M 2 peak is significantly shifted 
from its infinite value in Fig. [j]. At the same point (fc m ax) 
tends to zero in the infinite A limit, and tends to be a 
better indicator of the transition. All these facts are con- 
sistent with the interpretation of critical point 
of a phase transition characterized by the appearance of 
an "infinite" or percolating cluster. 



B. Disordered Systems 

Disordered systems are traditionally analyzed by con- 
sidering the distribution of an observable at a particular 
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point in parameter space. The important observables in be revealed by such a (potentially) blind test, 
disordered systems are Y = Y2 and fc max , and their distri- 
butions show significant non-Gaussian behavior, and this 
feature is in some sense a characteristic of the disorder. 

For the model given by Eq. (g), these can be obtained 
by a recursive procedure. Define x k — (p/) k -\/k\. Then 
the probability for a pattern separates into a term like 
Eq. ([3]) times a term which depends only on the multi- 
plicity. The partition function, Za — (x + j)a-i/A\, is 
then essentially a sum of terms contributed by each mul- 
tiplicity class, namely Za — J2 m x (m-i;i--y) Z A- Each of 
these terms can be computed for any Xk from the recur- 



sion raZ™ — J2k Xk ^A~k- Then the probability distri- 
bution n(fc max ) = AZa/Za as in the previous section. 

To determine the probability distribution U(Y), the 
partition function is broken into terms in which Y is 

fixed, i.e. Z A = Ey Z A> so that U ( Y ) = Z \I Z A- 
For computational purposes, it is convenient to con- 
tinue to break these partition functions into compo- 
nents where the multiplicity is fixed, so that Z\ = 



A. Nuclear Systems 

The corresponding analysis of fragmenting nuclei con- 
siders the distribution of the second moment Y and the 
size of the largest cluster fc max . Computing n(V) and 
n(/c ma x) is not complicated for models with fragmen- 
tation pattern given by Eq. (|J). We have already dis- 
cussed how to obtain II(fc max ) from the partition func- 
The distribution 11(F) is obtained 



IV A 



(jTi - 



-1,1-7)4/1 



Using x k as defined above, 



we can compute Z™' Y via the recursion mZ™' Y = 

k Xk^ A -k 

For many values of x and 7, these distributions are sig- 
nificantly non-Gaussian, for instance as seen in Fig. |^(a) 
which shows the distribution n(V). Firstly, one notes the 
large finite-size fluctuations seen in H(Y) when Y/A 2 > 
1/2. These disappear in the infinite limit. The pres- 
ence of cusps in IT(y) at Y/A 2 = 1,1/2 and 1/3 can 
be seen in figure |](a) and persists in the infinite limit. 
These cusps are not strictly an indicator of disorder, but 
the non-Gaussian nature of the distributions accentuates 
their presence. 

Figure ||(b) which shows n(fc max ) reveals that the dis- 
tribution of sizes of the largest cluster varies widely for 
different models. All these cases are likely to have sizable 
largest clusters, however there is a nontrivial distribution 
about fc max /A =1/2 for these models, due in part to the 
large number of patterns with fc max = A/2. The cusps 
at fc max I A — 1/2, 1/3 etc. arise from singularities asso- 
ciated with breaking the system into fragments of equal 
size. For example, below fc max /A = 1/2, there is no con- 
tribution from binary fragmentation. The sudden cutoff 
of this contribution leads to a cusp. 



V. CORRESPONDING ANALYSIS OF MODELS 

In this section we reverse our methods of analysis, 
treating nuclear models as if they were disordered sys- 
tems, and vice versa. Such an approach is warranted 
since both descriptions can be described without refer- 
ence to the fragmenting mechanisms, and therefore the 
methods should be insensitive to the source of the data 
and sensitive instead to the underlying property we are 
trying to isolate in the data (i.e. criticality and disor- 
der). Whether the methods are indeed selective should 



tion in section 
by breaking the contributions to the partition function 
Za into classes with fixed Y, Za — Yly zY a n( ^ a P" 
plying U(Y) 



As in section [VB, these 



AZ\ 



partition functions can be computed by recursion, i.e. 
- k kx k Zlzf. 
If we analyze the fragmentation distribution given by 
Eq. (H) with a typical Xk — x/k T with t > 2 by plotting 
the distribution on the largest cluster size and the second 
moment, we find over most of the range that the distribu- 
tions are essentially Gaussian centered about their expec- 
tation values. Figure |3| shows this. When x is small and 
the pattern is usually dominated by a single large frag- 
ment, the atypical events with smaller large fragments 
are distributed in a manner similar to a disordered sys- 
tem with a low value of x, only with a much diminished 
probability. This arises almost entirely from binary frag- 
mentation events. 



B. Disordered Systems 

The corresponding analysis of disordered systems ex- 
amines the expectation values of the reduced moments 
as a function of a changing parameter. We consider two 
cases given by Eq. (H) : the randomly broken interval with 
p(z) = x(l — z) x ~ l with x > (7 = 1), and the spin-glass 
model with 0<a:<l(a; + 7 = l). To facilitate com- 
puting the reduced moments for the models described by 
Eq. M) we have the following results for ensemble aver- 
ages of n k 



/ v = f A \ {x)A-k{l)k-l 

\kj (x + 7)^-1 
A' 

( nj (n k - S jk )) = x jlkl ( A _ j _ k y (Vi-l{l)k-i 
v Qc-7 + lU-j-fc 



(13) 



(x + i)A- 



(14) 



With these we can compute (Y s ) and its variance. Before 
we do so, it is interesting to note some properties of {n k ) 
itself. 

Since J2k = A, k(n k )/A is properly a probability 
distribution. In this case, the distribution has appeared 
before as a distribution for a replacement urn model in- 
troduced by Polya |39]. This correspondence was noted 



5 



by us before in the special case of 7 = 1 [Q. A central 
feature in that case was the appearance of a power law in 
the fragment distribution (n&) = 1/k for x = 1. Such a 
power law applies in this more general case, since in the 
large A limit, we can approximate (rife) by 



£(£+7) 1/ fc 
w ~ r(x)r( 7 ) a ^ a 



7-1 



(15) 



Setting x = 1, (n&) ~ 7^4 1 ~ 7 /c 7 ~ 2 which shows a power 
law falloff with exponent r = 2 — 7. 

If we now allow A — » 00, fc/v4 — » W, the dis- 
crete model becomes a continuum one. In this case, 
f(W) = lim^^oo A{tiwa) is the continuum distribu- 
tion of fragment sizes, with J c \Wf(W)dW = 1 and 

lim A ^ 00 (r s )/A 2 = Jl W s f(W)dW. Using Eq. (|l|), we 
arrive at 



r(x + 7) 
r(x)r( 7 ) 



(1 - w)*- 1 ^- 



(16) 



i.e., Wf (W) i s a beta distribution. Derrida and Flyvb- 
jerg pl] , |2^ , |33[ | , extending results on spin-glasses due to 
Mezard et al. Q] and other disordered systems, arrived 
at these distributions before. All the cases they consid- 
ered are in fact special cases of the above expression. 
Specifically, 7 = l,x = 1/2 reproduces the quenched 
random map, 7 = l,x > reproduces the randomly bro- 
ken interval and 7 = 1 — x reproduces the Sherrington- 
Kirkpatrick spin-glass model We have already seen 
why the first two models reduce to this model. The spin- 
glass result on the other hand defies a simple explanation. 
Distributions on the largest and second largest cluster 
also agree with the known spin-glass results, and so it 
appears that this model does indeed capture the cluster- 
ing properties of the spin-glass. 

Returning now to the question of moments, we see that 
in this case we can arrive at closed form solutions. En- 
semble averages for the multiplicity and its variance can 
be derived by applying Eqs . (|13|) and (|14|) and some com- 
binatorial identities (Eqs. (p31|) and (p32|)) to arrive at 



(to) = — I 1 

7-I 



(z + 7- 1) A J 



(17) 



Var(TO) = L , ( t ' i*U 



7 - 1 (x + j)a-i V 7 - 1 ( x + i)a-i 
1 (x-7 + 1)^ 



7-1 (* + 1)a-i J ' (18) 
Similarly, the mean and variance of Y can be determined 



(Y) = 



7 



I 1 X 

x + 7 \ A 7 



(19) 



2x 1 (A-l)(x + 1 + A)(x + 7 + A-l) 
1 ' A*(x + 1 ) 2 (x + 1 +l)(x + 1 + 2) • { ' 

From this the relative fluctuations in Y, f(Y) = 
(Y 2 )/(Y) 2 — 1 can be computed. They are important 



when determining whether a system is self-averaging [[L3| . 
Consider the case 7=1, where f(Y) w 2x / ((x+2)(x+3)) 
in the large A limit. For x ~ 1, the fluctuations are large 
and the system is said to lack self-averaging properties. 
For x oc A the fluctuations are small and f(Y) ~ 1/A as 
expected. This is true in general in this model, and for 
the spin glass case x + 7 = 1 leads to the known result 
(Y 2 ) = 1/3 (F)(l + 2(F)) II in the infinite A limit. 

The self-averaging property is not difficult to under- 
stand. At one extreme, the nucleus can be broken into 
individual nucleons so that Y/A 2 = 1/A; thus, Y varies 
inversely as the size of the system. At the other ex- 
treme, when the A nucleons are in one cluster Y/A 2 = 1. 
This behavior is analogous to the two extreme limits for 
the participation ratio in Anderson localization. For a 
localized state Y/A 2 ~ 1/A and for an extended state 
Y/A 2 <~ 1. A localized state corresponds to the local- 
ization of the electron wave functions around many in- 
dividual sites, while an extended state is the opposite 
case. A percolation description of localized versus ex- 
tended states can be found in Ziman ]l2| . Specifically, 
for site percolation an extended state includes an incipi- 
ent infinite cluster. 

Applying these results, we plot the various expecta- 
tion values in Fig. ^. From these graphs one can com- 
ment on the possibility of criticality in this model. It 
is well known that the x > 0, 7 = 1 case is not criti- 
cal in the traditional sense. The zeros of the partition 
function occur on the negative real x axis, well isolated 
from the relevant parameter domain. For the x + 7 = 1, 
< x < 1 model, the reduced moments appear to peak 
near x = 1/2, which suggests that x = 1/2 may be a crit- 
ical point for the model. However, Var(m) peaks closer 
to x = 0.9. At a true critical point one would expect this 
fluctuation to be maximal at the critical point. Since 
for any model the various expectation values must peak 
somewhere, the existence of a maximum is insufficient 
evidence for criticality. Further evidence is inferred from 
the fact that the partition function (x + j)a-i/A\ has 
zeroes on the negative real x + 7 axis in a complexified 
parameter space. These zeroes are isolated from .T+7 = 1 
real space, and suggests that a traditional critical point 
in which the zeros of the partition function encroach the 
thermodynamically relevant parameter domain is com- 
pletely absent in this model. 



VI. CONCLUSIONS AND OUTLOOK 

In this paper we have shown that the random break- 
ing of an interval process when suitably generalized en- 
compasses a wide range of physical models both in nu- 
clear fragmentation and in disordered systems. The fact 
that all these models are implicit to a particular process 
strongly recommends that the analysis of their properties 
should be the same. Traditionally however this has not 
been the case. The models have arisen separately, mo- 



G 



tivated by particular features which theorists wished to 
capture. In nuclear fragmentation, the feature to repro- 
duce is a specific criticality marked by the appearance 
of an "infinite" cluster. Interest has therefore focussed 
on reduced moments whose behavior changes with the 
appearance of such a cluster. In disordered systems the 
interest instead has been on whether such systems lack 
self-averaging properties, i.e. that have significant non- 
Gaussian fluctuations. In this case, the focus is on prob- 
ability distributions of the measures which might express 
such non-Gaussian fluctuations, e.g. the largest cluster 
and the second moment. 

Given then this separate emphasis, it is understood 
why the same variables in the models are studied in dif- 
ferent ways. The corresponding analysis taken here rein- 
forces the distinct approaches, as it shows that the gen- 
eral behavior of the two classes of models differ consider- 
ably. Nevertheless, the convergence of their description 
within the random breaking of an interval mechanism 
makes such a distinction somewhat puzzling. This mech- 
anism is unlikely to be physical in the nuclear case, but 
this is not an explanation for the different behavior in 
the two cases. Indeed, the explanation may have a great 
deal to do with criticality or its absence. Models with ex- 
ponent t < 2 tend to look disordered, r > 2 tend to have 
distinct critical points and Gaussian distributions. Since 
t = 2 + 1/(5, the two regions are characterized by positive 
and negative 5. It would be interesting to confirm this 
observation as a general feature. 
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the fc objects and (n — fc)! ways of permuting the re- 
maining elements which are not in this cycle. So the 
probability is (fc - l)!(n - fc)!(^I 1 1 )/n! = 1/n. What 
about the second cluster? The (n — fc)! ways of permut- 
ing the elements not in the first cluster is a permuta- 
tion p' G 5* n _Ai chosen uniformly at random. In this 
permutation, A2 is the first element, and by the above 
argument has size X2 — k with probability l/(n — Ai). 
Iterating yields the probability of the whole sequence as 
being Pr(Ai, . . . , A m ) = l/(n(n - Ai)(n - Ai - A 2 ) • ■ •)■ 
Summing this over all possible orders one obtains l/n! 
times the number of permutations with this particular 
cycle class structure. It is well known that this is given 
by Cauchy's formula [Q n!/ J7 fc riklk nk , so that 



E 



1 



TT-f 1 

A , A n k \ \k 



n(n — Ai)(n — Ai — A2) ■ 



(AI) 



APPENDIX B: SUM RULES 

For systems satisfying Eqs. (|l^) and (|l4|) one has the 
following sum rules 



k—p 



[A] 



p+q 



(rij(nk - Sjk)) = x 



(Bl) 

ph^lph^l , (B2) 

{x + Vp+q-l 



where p, q > 0, [fc]o = 1, [k] p = k(k — 1) . . . (fc — p + 
1). Notice that the right hand sides are independent of 
A. These rules are helpful in determining expectation 
values of m, Y etc. and can be derived using a common 
combinatorial identity known as Norlund's formula p3|. 



APPENDIX A: A COMBINATORIAL IDENTITY 



Recall that a permutation on n elements can be de- 
scribed by the cycles containing the elements which 
permute among themselves under the iterated action 
of the permutation. For example, the permutation 
P = (3125468?) * s usuau y written in cycle notation as 
(132) (45) (6) (78). If each cycle of length fc is associated 
with a cluster of the size fc, and these clusters are or- 
dered by the largest element in each cluster, then each 
permutation can be mapped to a sequential partitioning. 
The above example would correspond to the sequential 
partition (2, 1, 2, 3). 

Suppose one chooses a permutation p 6 S n uniformly 
at random, i.e. with probability l/n!. What is the 
probability that the first cluster in the sequence has size 
Ai = fc? There are (tZi) ways of choosing the fc — 1 
elements which are not element n but are in n's cycle 
map. There are (fc — 1)! different cycles which permute 
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FIG. 1. The expected value (a) and variance (b) of the multiplicity and the expected values of the reduced second moment 
(c) and the largest cluster (d) for the model given by Eq. (H) with Xk = x/k T , r = 2.5 (solid line), r = 3.0 (dashed line). 
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FIG. 3. The probability distribution of (a) the second moment 11(F) and (b) the size of the largest cluster Il(A: max ) for the 
model given by Eq. (|) with x k = x/k 5/2 for various x. 
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